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^ ■ We present results from numerical studies of the NJL model with two massless quarks 
at nonzero temperature. We show that the model undergoes a second order chiral phase 
^ ! transition which belongs to the 3d 0(4) spin model universality class. 
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1. INTRODUCTION 

Massless QCD with two flavors undergoes a chiral symmetry restoration transition 
which is believed to be of second order. According to arguments by Pisarski and Wilczek 
[H] the universality class of the transition is that of the dimensionally reduced 3d 0(4) 
Heisenberg spin model. The reasoning behind their conjecture is that quantum fluctu- 
ations become irrelevant near the critical temperature T c . As T — > T c the correlation 
length £ diverges and the characteristic mode is uo c ~ 0. The fermions are expected to 
decouple from the scalar sector because they do not have zero modes. The energy of the 
zero point quantum fluctuations of characteristic frequency lu c , Tiuj c is much less than the 
thermal energy per degree of freedom fcgT. Therefore, in this regime one observes the 
critical behaviour of a classical (thermal) transition. However, there is a possible loophole 
in this scenario. In fermionic models the mesons are composite // states whose density 
and size increase as T — > T c . This may imply that the fermionic sub-structure is apparent 
on physical length scales and as a consequence we may have a maximal violation of the 
bosonic character of mesons near T c . If this is the case then the quarks become essential 
degrees of freedom irrespective of how heavy they are. 

Since the problem of chiral symmetry breaking and its restoration is intrinsically non- 
perturbative the number of techniques available is limited and most of our knowledge 
about this phenomenon comes from lattice simulations. Various lattice results indicate 
that the QCD chiral phase transition at T ^ is second order However, due to 

problems with finite size effects near the critical temperature and the proper inclusion of 
fermions on the lattice, these calculations are still not fully satisfactory. Therefore, an 
important role is played by model field theories which incorporate various basic features 
of QCD and they are much easier to deal with. Such a theory is the Nambu— Jona-Lasinio 
(NJL) model. The NJL model was introduced in the sixties as a theory of interacting 
nucleons and later it was reformulated in terms of quark degrees of freedom. In the 
continuum, the model is described by the Euclidean Lagrangian density 
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where r = (ri, T3) are the Pauli spin matrices which run over the internal SU(2) isospin 
symmetry. The index i runs over N c fermion colours and g 2 is the coupling constant of 
the four-fermion interaction. The model is chirally symmetric under SU(2)l ® SU(2)r\ 
i/j — > (PlU + PjtV)ip, where U and V are independent global SU{2) rotations and the 
operators Pl,r = |(1 ± 75) project onto left and right handed spinors respectively. It 
is also invariant under U(l)v corresponding to a conserved baryon number. The theory 
becomes easier to treat, both analytically and numerically, if we introduce scalar and 
pseudo-scalar fields denoted by a and 7? respectively. The bosonised Lagrangian is 

C = ^ (0 + a + • r) ^ + ^ (V + vr • vr) . (2) 

The dynamic generation of fermion masses brought about by the breaking of chiral sym- 
metry at g 2 larger than a critical g 2 is one of the important properties of the model. This 
implies the creation of the Nambu-Goldstone bosons (n) and the creation of the chiral 
partner sigma (cr) which are // bound states. The interaction strength has a mass dimen- 
sion -2, implying that the model is non-renormalizable. It works well in the intermediate 
scale region as an effective theory of QCD with a momentum cutoff comparable to the 
scale of chiral symmetry breaking A. x $b — lGeV. Various authors demonstarted the suc- 
cess of the model as a quantitative theory of hadrons (for a review see e.g. [ 4J). It was 
also shown that the model has an interesting phase diagram in the (temperature, quark 
chemical potential n) plane, in close agreement with predictions from other effective field 
theories of QCD. In the Hartree-Fock approximation, the theory has a tricritical point 
in the (T, fi) plane [ E] when the quark current mass, the coupling g and the momentum 
cut-off A are chosen to have physically meaningful values. Recently, lattice simulations 
showed that the ground-state at high /1 and low T is that of a traditional BCS superfluid [ 
|U] . It was also shown numerically that the d- dimensional Z 2 symmetric model undergoes a 
finite temperature phase transition which belongs to the (d — l)-dimensional Ising model 
[ 13 IE] implying that the dimensional reduction scenario is the valid description of the 
thermal transition. Here we study the case where the chiral symmetry is continuum. 

2. RESULTS 

We performed lattice simulations of the NJL model with the Hybrid Monte Carlo 
method. Details about the staggered fermion action used in our simulations can be found 
in [E]. The asymmetric lattice has L s lattice spacings a in spatial directions, L t <C L s 
lattice spacings in the temporal direction and volume V = L\L t . The temperature is given 
by T = l/(L t a). On a finite volume and with the quark bare mass set to zero, the direction 
of symmetry breaking changes over the course of the run so the chiral condensate averages 
to zero over the ensemble. Another option is to introduce an effective order parameter 
$ = y\J(J2 x c(x)) 2 + J2i(J2x 7r i( a; )) 2 ) which in the thermodynamic limit is equal to the 
true order parameter £ = (a) extrapolated to zero quark mass. In order to study the 
critical behaviour on lattices available to us we use the finite size scaling (FSS) method, 
which is a well-established tool for studying critical properties of phase transitions. On a 
finite lattice the correlation lenght £ is limited by the size of the system and consequently 
no true criticality can be observed. The dependence of a given thermodynamic observable, 
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P, on the size of the box is singular and according to the FSS hypothesis is given by: 

P(t,L s )=L^Q P (tLl/n, (3) 

where t = (/3 C — f3)/[3 c , v is the standard exponent of the correlation length and Qp 
is a scaling function, which is not singular at zero argument. The exponent pp is the 
standard critical exponent for the quantity P. Studying the dependence of P on L s and 
using eq. (j3J) one can determine such exponents. We performed simulations with N c = 8, 
fixed L t = 4 and L s = 8, 54 at three or four values of (3 for each lattice size and then 
used histogram reweighting techniques to calculate the values of the observables at nearby 
couplings. The temperature T is varied by varying the lattice spacing a(g). Near (3 C the 
expansion of the FSS scaling relation for the susceptibility x — V(§ 2 ) is: 

X = (do + ait + alt 2 + ...) + (6q + b x tL 1 J v + b 2 t 2 L 2 J v + ...)£ 7/ ". (4) 

For large enough L s one can neglect the terms in the first parenthesis. By fitting our 
data to eq. Q we get: (3 C = 0.52579(3), v = 0.74(3), j/u = 1.95(1). These results are in 
agreement with the 3d 0(4) spin model exponents [ 9J v = 0.748(9) and 7/V = 1.9746(38). 
A method usually used to determine (3 C is by measuring the Binder cumulant, 

As a consequence of eq. (JHI), near /3 C , Ul = f(tL l l u ). In the large L s limit Ul has a 
universal value U* at (3 C . Close to p c , U(L s ,f3) ~U* + a x {l3 c - (3)L l ' v + a 2 {(3 c - (3) 2 L 2 / V . 
The value of U* in the classical 3d 0(4) universality is 0.636(1). Fitting the above relation 
to data generated on lattices with L s = 30,40,54 we get v = 0.69(4), U* = 0.616(2) and 
(3 C = 0.5261(1). The measured values of v and U* are close to their respective values in 
the 3d 0(4) spin model. 

It has been shown analytically that in the infinite iV c limit the theory undergoes a 
second order phase transition with Landau- Ginzburg mean field (MF) exponents for any 
dimensionality d [EH]- The reasoning behind this behaviour is that bosonic fluctuations 
are neglected, whereas the fermions decouple because their non-zero Matsubara frequency 
acts as a mass term of order irT for the (d — l)-dimensional components of the field. 
However, as we showed in the previous paragraph for iV c = 8 the dimensional reduction 
scenario is the correct description of the transition. It has been shown in [|7j that although 
the fermions do not affect the actual critical behaviour they do affect the way the universal 
behaviour sets in. Close to T c , MF theory is expected to break down and in the case of 
d — 4 the non-trivial scaling region is suppressed by a factor ~ N~ x [|7j. 

We perfomed simulations in the broken phase on 4 x 54 3 lattices with N c = 8 and 
N c = 16. In the MF region the order parameter E is equal to the fermion mass and 
therefore acts as an inverse correlation length. In Fig. ^ we plot S 2 versus (3. The data in 
the MF region fits well to straight lines because in this region the critical exponent for S 
is /3 mag = 1/2. For N c = 8 the quality of the linear fit deteriorates significantly when we 
include the (3 = 0.520 point. For N c = 16 all the (3 > 0.515 data give a high quality fit, 
implying that in this case the crossover from the MF region to the 3d 0(4) region occurs 
at larger correlation lengths than for N c = 8. This result verifies that the 3d 0(4) scaling 
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Figure 1. S 2 



vs (3 for N c = 8, 16. 



Figure 2. £ vs (3 close to (3 C for N c = 8. 



region is suppressed at large N c . We also fitted the N c = 8 data for (3 = 0.520 — 0.5225 
(see Fig. |2J) to the scaling relation S = const. (/3 C — fl) 13 ™^ with fixed /3 C = 0.5260 and 
extracted /? mag = 0.360(15) which is close to the 3d 0(4) exponent /3 mag = 0.3836(46) 

3. SUMMARY 

We showed numerically that the NJL model of QCD with two massless quarks un- 
dergoes a second order chiral phase transition which belongs to the 3d 0(4) spin model 
universality. Furthemore, as shown some time ago in the (2 + l)d Z 2 -symmetric Gross- 
Neveu model [|7j the actual critical region is suppressed if the number of fermion species 
in the system is large. It will be interesting to check with comparable precision whether 
this scenario is the correct description of the QCD thermal phase transition. 

Discussions with S. Chandrasekharan, S. Hands and J. Kogut are greatly appreciated. 

REFERENCES 

1. R. Pisarski and F. Wilczek, Phys. Rev. D29 (1984) 338. 

2. F. Karsch, Phys. Rev. D49 (1994) 3791; F. Karsch and E. Laermann, Phys. Rev. 
D50 (19 94) 6954; S. Aoki et al, Phys. Rev. D63 (2001) 034502; J.B. Kogut and D.K. 
Sinclair, |arXiv:hep-lat/0408003| 

3. Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 345; Phys. Rev. 124 (1961) 
246. 

4. T. Hatsuda and T. Kunihiro, Phys. Rep. 247 (1994) 221. 

5. T. Hatsuda and T. Kunihiro, Phys. Lett. B198 (1987) 126. 

6. S.J. Hands and D.N. Walters, Phys. Rev. D69 (2004) 076011. 

7. J.B. Kogut, M.A. Stephanov and C.G. Strouthos, Phys. Rev. D58 (1998) 096001. 

8. S. Chandrasekharan and J.C. Osborn, Phys. Lett B496 (2000) 122; S. Chandrasekha- 
ran, J. Cox, K. Holland and U.J. Wiese, Nucl. Phys. B576 (2000) 481; J. Cox and K. 
Holland, Nucl. Phys. B583 (2000) 331. 

9. K. Kanaya and S. Kaya, Phys. Rev. D51 (1995) 2404. 

10. A. Kocic and J.B. Kogut, Phys. Rev. Lett., 74 (1995) 3109. 



